Matter-wave interference in s-wave and p-wave Fermi condensates 
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We discuss the time evolution and matter-wave interference of Fermi condensates on the BEC 
side of Feshbach resonances for s and p-wave superfluids, upon release from harmonic traps. In s- 
wave systems, where the order parameter is a complex scalar, we find that the interference patterns 
depend on the relative phase of the order parameters of the condensates. In p-wave systems involving 
the mixture of two-hyperfine states, we show that the interference pattern exhibits a polarization 
effect depending on the relative orientation of the two vector order parameters. Lastly, we also 
point out that p-wave Fermi condensates exhibit an anisotropic expansion, reflecting the spatial 
anisotropy of the underlying interaction between fermions and the orbital nature of the vector order 
parameter. Potential applications of our results include systems of ultra-cold atoms that exhibit 
p-wave Feshbach resonances such as 6 Li or 40 K. 

PACS numbers: 03.75.Ss, 03.75.-b, 05.30.Fk 
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I. INTRODUCTION 

Matter-wave interference is a very powerful tool to 
study quantum phase coherence between atomic Bose 
Einstein condensates (BEC) [ESQ, and spatial quan- 
tum noise of bosons in optical lattices [4j. Similar 
techniques can also be applied to study Fermi conden- 
sates 0, H, 0, H, H, [HI , where superfluidity can be tuned 
from the BCS to the BEC regime. These experiments 
may reveal that the time dynamics in the BCS regime 
is overdamped (large Cooper pairs can decay into two 
atoms), while in the BEC regime it is essentially un- 
damped (tightly bound molecules are stable) [ll|, fl2| |. 
Matter-wave interference experiments of s-wave Fermi 
condensates may be readily performed, since stable con- 
densates already exist. For s-wave Fermi condensates 
in the BEC regime quantum interference effects are ex- 
pected to be similar to those of atomic Bose condensates, 
and the interference pattern should depend essentially on 
the phase difference of the order parameters between two 
interfering clouds. 

In contrast, it is more interesting to study interference 
effects in p-wave superfluids because of the vector nature 
of the order parameter. Many groups have reported some 
progress towards the formation of p-wave Fermi conden- 
sates in harmonically trapped clouds [H, [lj, [HI, [l6| and 
in optical lattices [17| , where p-wave Feshbach resonances 
have been observed. Similar to the s-wave case, the 
scattering cross section for p-wave collisions has a peak 
that rises over three orders of magnitude above the small 
background cross section near p-wave resonances, thus 
suggesting a divergence of the scattering parameter (vol- 
ume). Therefore, matter- wave interference can be po- 
tentially observed using Feshbach resonance techniques 
when two-body dipolar or three-body losses are not too 
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large (l7j . 

In addition to the diverging scattering parameter (vol- 
ume), p-wave Feshbach resonances also reveal character- 
istic features in contrast to s-wave counterparts. These 
features include splitting of resonance peaks depend- 
ing on hyperfine (pseudospin) states as in 6 Li and 40 K 
(i.e., |11), 1 12) + |21), and |22)) [H, or 

splitting depending on angular momentum projections 
as in 40 K ( i.e., the magnetic quantum number mg — 
or ±1) [lj, [13. These experimentally observed split- 
ting allow the possibility of tuning separately p-wave 
scattering parameters in different pseudospin and/or mi 
states, such that the p-wave interaction can be in general 
anisotropic in both pseudospin and angular momentum 
states. These additional degrees of freedom are respon- 
sible for the much richer matter wave interference phe- 
nomena in p-wave Fermi condensates than in the corre- 
sponding s-wave case. 

Prior and in parallel to experiments the BCS to BEC 
evolution of p-wave superfluids was discussed in the con- 
text of ultra-cold atoms 0, GJ, [H HH, [H, H, H , where 
the existence of quantum phase transitions was empha- 
sized OESIMliilSi]. However, these works dealt only 
with zero and finite temperature thermodynamic proper- 
ties. 

In this manuscript, we discuss the time evolution and 
matter-wave interference of s-wave and p-wave Fermi 
condensates on the BEC side of Feshbach resonances, and 
extend our previous work in this area [25] . Our main re- 
sults are as follows. While in atomic BEC and s-wave 
Fermi superfluids quantum interference patterns depend 
essentially on the relative phase of the two clouds, we 
find that in p-wave Fermi superfluids there can also be a 
strong dependence on the relative angle between the two 
vector order parameters, thus producing a polarization 
effect. This polarization effect is a direct consequence 
of macroscopic quantum coherence of a large number of 
molecules, as well as the vector nature of the order pa- 
rameter. We also discuss the Josephson effect between 
two p-wave condensates and show that it depends not 
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only on the existence of phase coherence, but also on the 
relative orientation of the vector order parameter of the 
two clouds. Furthermore, we also show that p-wave Fermi 
condensates exhibit an anisotropic expansion, reflecting 
the spatial anisotropy of the underlying interaction be- 
tween fcrmions and the orbital nature of the vector order 
parameter. 

The remainder of this manuscript is organized as fol- 
lows. In Section HH we discuss the dynamics of Fermi 
condensates with s and p-wave interactions. In partic- 
ular, we derive the equation of motion for the vector 
order parameter of p-wave condensates in the strongly 
interacting BEC limit. This equation of motion is fur- 
ther discussed in Section ITTTl where the time-of-fllght ex- 
pansion of a harmonically trapped cloud is studied. In 
Section llV| we describe matter- wave interference of two 
condensates, and demonstrate that the interference pat- 
tern for the p-wave case depends crucially on the relative 
orientation of the two vector order parameters. In Sec- 
tion [V] we relax the restriction of being near the BEC 
limit by moving towards unitarity and demonstrate that 
the time-of-fiight expansion of a p-wave Fermi conden- 
sate is in general anisotropic. Lastly we summarize our 
main results in section IVII 



II. 



EFFECTIVE THEORY AND EQUATION OF 
MOTION 



We consider a system of fermions with mass m in two 
hyperfine states (pseudospins), labeled by greek indices 
a = 1,2. The Hamiltonian density is (with h = ks = 1) 



V 2 

~ + I7«t(r, t) 
2m 



■0a (r, t) 



W(r,*)=V4(r,i) 



dr' 



where repeated greek indices indicate summation, tpj^ 
{ipa) are creation (annihilation) operators of fermions in 
state a, t/ ex t(r,t) is the time dependent trapping poten- 
tial, and V a /3-y5 = V Q /3 7 5(r — r'). The Hamiltonian of the 
system is then H(t) = J drH(r, t). 

The generating functional for non-equilibrium pro- 
cesses associated with H(t) is [26| 

Z(t) = Trrt(t,t )exp[-f3(H(t ) - n a N a )}U(t,t ), (2) 

where N a is the number operator for fermions of type 
a, \i a is the corresponding chemical potential, and 
f3 = l/T is the inverse temperature. Here, U(t,to) = 
cxp[— i f tg H(t')dt'} is the time evolution operator. This 
expression implicitly implies that the trapping potential 
C/ ox t is time independent for t < 0, such that the ini- 
tial condition corresponds to a thermal equilibrium state. 
This assumption is directly related to experiments where 
the traps are effectively static before the cloud release. 
Therefore, the generating functional at any time to < 
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FIG. 1: Integration contour C used in Eq. Q. 



takes the form Z(t ) = Trexp [-0(H(t o ) - fi a N a )}, 
which is just the partition function describing thermal 
equilibrium properties. 

By introducing a complex time r, the generating func- 
tional Eq. @ can be written as 

Z(t)= [ V[i)l,^ a }ei- S2 ^ a] - Si[ ^'^}, (3) 
Jbc 

where the boundary condition (BC) is antisymmetric 
over the integration contour C of the functional integral 
(as shown in Fig. [1]), and S2 and £4 are quadratic and 
quartic action functions of fermionic operators ip, respec- 
tively. 

In what follows, we discuss the p-wave case in detail 
and quote the more standard results for s-wave. To 
calculate the effective action of the system, we write 
Va/3-ys{p) = V(p)T a p 1 s, where in a triplet channel 
r a /37(5 = v Q( 3 • (v')yg. The pseudospin matrix elements 
for the j th component of v are (vj) a /3 = (io'j&y)a0, where 
<7j are Pauli matrices. By writing the interaction in this 
form, we implicitly assume that the interaction is sym- 
metric in pseudospin space with Vim = V1212 = V2222- 
It should be emphasized that this assumption is intro- 
duced just for simplicity and the resulting equation of 
motion is valid for a more general interaction. In fact, 
the experimental realizations of the polarization effect in 
mater-wave interference and the anisotropic expansion 
proposed here rely on the tunability of the interaction 
in both pseudospin and angular momentum (rrii) spaces. 
Here, we first consider the simpler case of symmetric in- 
teraction, and postpone the general discussion until later 
in Section [TV] 

In order to integrate out the fermions, we introduce 
the field B^(r, r',r) = Vi(r, r)v a/ 30^(r', t) and the cor- 
responding auxiliary field d(r,r',T). Since the d field 
depends on two spatial variables r and r', it can also be 
transformed into the center-of-mass R = (r + r')/2 and 
relative p = r r' coordinates: 



d(r,r',r) = 



E 



D 



(R,r)r?„ Am ,(p), (4) 



where rj n ,£ iTne (p) are eigenfunctions of the reduced two- 
body Hamiltonian H2 — — V 2 /m + V(p). 
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The vector nature of the order parameter d or D can 
be understood from the pseudospin structure of the pair 



wave function ty z 



ffn|ll)+ffi 2 (|12) + |21)) + ff22 |22) 



Using the symmetric matrices iaja y , this pair wave func- 
tion can be represented in the form ^ 
where the components of d are related to the amplitudes 



pair 



> y ■ 



through .gn 



-idy, 512 = d z , and 322 = d x +id y . 



Therefore, the direction of d (or D) determines the am- 
plitude of the pair wavefunction in each of the pseudospin 
triplet channels |11), 1 12) + |21), and |22). For 6 Li these 
orthogonal states correspond to the v — 38 vibrational 
state of the singlet potential [l6[ of total electronic spin 
S = 0, and total nuclear spin 1=1. 

For dcfinitcness, we consider a pure p-wave interac- 
tion where the ground state is three-fold degenerate 
(£ = 1) and labeled by me. = —1,0,1. A rotation 
of basis from spherical harmonics Yi >nu {p) to p v =x.y,z 
with corresponding cigcnfunctions »7o,i(p)Pi/ leads to d = 
Y^, v T )Vo,i(p)Pu a t l° w temperatures, where the 

higher energy states are not excited. In the BEC limit, 
Cooper pairs are tightly bound molecules and the rela- 
tive degrees of freedom p can be integrated out leading 
to an effective action 



.S',„ = -J dr J dB.{Dt(R,T)- £D(R,r) 
[2|D(R,t)| 4 -|D 2 (R,t)| 2 ]}, 



(5) 



where -Dj(R) = H„ -D;>(R), and the operator K. = 
id T — 2{/ cx t(R, t) + V^/ (4m) corresponds to the action of 
an ideal non-equilibrium gas of Bose particles with mass 
M = 2m. This action leads to equations of motion 



id t D 3 = 



-J| + 2lWR,i) + 2 3o |D| 



D, 



-g (T>-T>)D]. 



(0) 



Notice that this expression is different from the time- 
dependent Gross-Pitaevskii (TDGP) equation for an 
atomic vector boson held. The difference comes from the 
last term, which describes a non-unitary complex order 
parameter of the underlying paired fermions. In contrast, 
the standard TDGP equation for scalar atomic bosons is 
obtained in the s-wave case [26| . 

Equation © can be simplified to the TDGP form in 
two special cases. First, if the atomic hyperfine states 
|1) and 1 2) are equally populated with N\ — N2 = N 
{fix = (X2 = A*) arid D is unitary, then D is a real vector 
with an overall phase, leading to the equation of motion 



id t Dj 



2M 



■2U ext (R,t)+g \T>\ 



D, 



(7) 



Second, if only one atomic hyperfine state is populated, 
then D is non-unitary and D = A(l,±i,0), where A 
is a complex constant. Thus, the last term D • D in 
Eq. ^ vanishes, and the equation of motion is identical 
to Eq. 0, with go — > 2go. In the following section, we 



confine ourselves to these two special cases and discuss 
the time-of-flight expansion of a Fermi condensate upon 
released from the trap. 



III. TIME-OF-FLIGHT EXPANSION OF A 
HARMONICALLY TRAPPED CLOUD 

In order to describe the time evolution of a triplet 
Fermi condensate, one needs to solve the nonlinear equa- 
tion of motion Eq. ([7|), which usually requires special 
numerical treatment. For simplicity and definiteness, we 
consider a Fermi superfluid is released from a harmonic 
trap, i.e., 



U cxt {r,t)= ma#)rf/2, 



(8) 



3=x,y,z 



where u>j(t < 0) = u>j are constants and u>j(t > 0) 
Thus, for t < 0, the system is described by 



Mo-Dj(R) 



-|^ + 2C/ cxt (R)+. g |D| 2 



Dj, (9) 



where /io is the effective boson chemical potential, and 
.9 = 9o (9 = 2<7o) when D is unitary (non- unitary). For 
dominant Boson interactions the Thomas-Fermi approx- 
imation leads to 



|D(R,0)| 



^o-2U^ (R) i for Mq > 2 [/ cxt (R) (1Q) 
0, otherwise. 



When this approximation fails the initial condition for 
the time evolution can be obtained by solving Eq. ([5]) 
numerically. 

For t > 0, we use the transformation RAt) = 
bj(t)Rj(0), where the scaling factors bj(t) satisfy [13, HI], 



dt 2 



A(t)bj(t) 



(11) 



with A(t) = b x (t)b y (t)b z (t) and initial conditions bj(0) = 
1. For a cigar-shaped trapping potential with axial sym- 
metry (u> x = oj y = u>j_ 3> u> z ), Dj (R, t) becomes 

^t) M ^ ^0), (12) 
vl + * 

where Rk = Rk/bkif) are scaled coordinates, A = uj±t is 
the dimcnsionlcss time, and the phase factor 



Rl(t) db k (t) 
b k (t) dt 



(13) 



The result for s-wave is formally identical to that of 
Eq. (I12p with the substitution Dj — > 1 J r , where \t rep- 
resents the scalar order parameter. 



In the limit where e = ui z /uj± <C 1, approximated so- 
lutions for bj (t) can be obtained as a power expansion of 
e, leading to 



S(R(t),t) -- 



Ho tan 1 (X) 

LU I A 



0(A) 



+M 



1 + A 2 



{R 2 X + R 2 V ) + Me 2 Lu A _tan- 1 (X)R 2 zl (14) 



where the function 
following form 



9(A) in the second term takes the 



6(\) 



Notice that 8(\) 
term in Eq. (fT4")l 

when tln(A) ~ e 
ters where e ~ 10 



2a; tan" 1 (sc) - ln(l + x 2 ) 



2(1 

ln(A) as t 



+ x 2 ) 



dx. 



(15) 



oo, hence the second 
becomes comparable to the first one 
For realistic experimental parame- 
2 and ui± ~ 10 3 Hz, this condition is 



satisfied only when t ~ 10 s, which is several orders of 
magnitude larger than any time scales in current exper- 
iments 0, 1!]. Therefore, the 0(A) term in Eq. fH} can 
be neglected in the following discussion. 

It should be emphasized that the scale transformation 
discussed above is valid under the hydrodynamic approx- 
imation, where the collisions between molecules are not 
playing a crucial role. Furthermore, it is also assumed 
that the Cooper pairs (Feshbach molecules) are long lived 
and the damping processes are negligible. With these 
conditions in mind, we discuss next matter wave inter- 
ference of two clouds. 
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FIG. 2: Interference pattern versus dimensionless time A = 
oj±_t for p-wave Fermi condensates in the BEC limit with 
lj z /lo± = 0.1, assuming |D£ ■ Dr| is maximal. The plots 
include only the superfluid part, and show columnar density 
versus x, y coordinates in units of the initial clouds separation 
W. The patterns are similar to those of atomic scalar bosons, 
and s-wave paired Fermions. 



IV. MATTER- WAVE INTERFERENCE AND 
POLARIZATION EFFECT 

In this section, we consider first the matter- wave inter- 
ference of two spatially separated condensates such that 
the energy barrier between them is large enough to ne- 
glect the tunneling effect. Thus, one may write down the 
total wave function as 

$ tot (R,<) = $ L (R,f) + $ R (R,i), (16) 

where $p cc i ■ D^-pOjO y denotes the pair wavefunction 
of Fermi condensate in the left (P = L) or right (P = R) 
trap. The left and right trap centers lie at (— W/2, 0,0) 
and (W/2, 0,0), respectively, where W is the distance 
between traps. We consider the case of two identical 
axially symmetric Fermi condensates with to± 3> u) z . In 
this case, the time evolution of each cloud is described 
by 



D j>P (R,t) 



exp[iS(R±W±/2,t)] 



D J - )P (R± Wx/2,0), 



VT+A^ 

Thus, for a single run of experiment, the particle density 
n(R,t) = |$tot(R,i)| 2 is 

n(R,t) cx |D L (R,i)| 2 + |D R (R,i)| 2 

D[(R + TFx/2, 0) • Dp, (R - VFx/2, 0)e« 



-2Re- 



A(X) 



-,(17) 



where the time dependent phase factor x(R-) i) — S(R + 
Wx/2,t) - S(R - Wx/2,t) + xo, and X o is the initial 
relative phase of the two condensates. The result for s- 
wave is formally identical to that of Eq. (TPTj) with the 
substitution Dj — * ty. 

When each cloud has the same hyperfme state occupied 
(e.g. ^Ppair = <?ii|ll)) the D vectors in each cloud have 
the fixed form A(l,i,0) and fringes are present in all 
experimental realizations. This result is similar to the s- 
wave case where the order parameter is a complex scalar. 
However, when both Fermi condensates are in unitary 
states, D is essentially a real vector with an overall phase, 
and n(R, t) shows an angular dependence controlled by 
the dot product term in Eq. (fl~7)) . When the two order 
parameters are parallel, this term is maximal and the 
interference pattern is most visible (Fig. [5]). But if the D 
vectors are perpendicular, fringes are absent at all times 
(Fig. [3]). Therefore, in the unitary case the existence and 
intensity of interference fringes are very sensitive to the 
relative orientation of the D vectors. 

This sensitivity to the relative orientation of the or- 
der parameters and the corresponding polarization ef- 
fect also manifest themselves in the Josephson tunneling 
between two condensates. By considering a tunneling 
process across the energy barrier, the left and right con- 
densates can be described by the modified equations of 
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FIG. 3: Interference pattern versus dimensionless time A = 
u>±t for p-wave Fermi condensates in the BEC limit with 
uj z /uj± = 0.1, assuming |D[ ■ D R | = 0. The plots include 
only the superfluid part, and show columnar density versus 
x, y coordinates in units of the initial clouds separation W. 



motion 
id t D 



-|| + 2C7 ext (R)+ ff |D L | 



J tffl. / r J y(R,R / )DRj(R / ,t), (18a) 



■|^ + 2C/ oxt (R)+ 5 |D R | 2 



D 



R,j 



I dR'T* k (R,R')D LJ (R',t), (18b) 



where X^R, R') = T c + T s is the tunneling matrix, 
with T c and T s are the pseudospin-conserving and non- 
conserving portions, respectively. In the case where the 
trapping potential is much smaller than the energy dif- 
ference between the two hyperfine states, or the time 
scale of experiments is much shorter than the lifetime of 
each hyperfine state, the non-conserving tunneling pro- 
cess is negligible, leading to Tjk — T(R, R')Sjk- Thus, 
the Josephson current J is dominated by the pseudo-spin- 
conserving tunneling processes, leading to 

J = Im / dR [ dR"T(R,R')D[(R,i) . D R (R',i). 
Jls J 

(19) 

Notice that the Josephson current depends crucially on 
the relative orientation of the two vector order parame- 
ters, and acquires a polarization effect as in matter-wave 
interference. In particular, if D^(R) ■ D R (R') = is sat- 
isfied, there is no Josephson tunneling between the two 
condensates for all R and R'. This property is similar to 



that encountered in Josephson junctions of lattice p-wave 
superconductors, where the Josephson tunneling current 
between two p-wave crystals depend on the relative ori- 
entation of the vector order parameters [29|, [3(| • 

Since the interference pattern and the Josephson cur- 
rent depend crucially on Dl and D R , it is important 
to understand how these vectors can be controlled ex- 
perimentally. In the discussion above we assumed a 
symmetric interaction in pseudospin space, i.e., Vim, 
V^i2i2 and V2222 were identical. However, experimental 
results for p-wave Feshbach resonances show a finite sep- 
aration in different channels [16| implying different in- 
teractions in pseudospin space. When the different in- 
teraction strengths are absorbed into an effective D vec- 
tor, an equivalent procedure leads to equations similar 
to Eqs. ^ and ((7|). For instance, the p-wave reso- 
nances for 6 Li occur at 159G (width 0.4G), 185G (width 
0.2G) and 215G (width 0.4G) for the |11), 1 12) + |21), 
and 1 22) channels, respectively. By applying a constant 
plus a gradient magnetic field, the local field at L (R) 
cloud can be tuned to be 216 (214)G, which is above 
(below) the 1 22) resonance. Thus, a sweep down of the 
constant magnetic field by 30G makes the L cloud cross 
the 1 22) but not the 1 12) resonance, while it makes the 
R cloud cross only the 1 12) resonance. In this case, the 
L cloud is in the BEC regime of the 1 22) channel, with 
* P air,L » 522,l|22) or D L = 2 2,l(1/2, -i/2, 0). However, 
the R cloud is in the BEC regime of the 1 12) channel, with 
IWr « <?12,r(|12) + |21» or D R = <?12,r(0, 0, 1). After 
these initial states are prepared the total magnetic field 
B(r) (constant + gradient) and the harmonic trapping 
potential are turned off suddenly (very fast) and simul- 
taneously, as it is standard in experiments pH |32| , such 
that the magnetic field does not interfere with the sub- 
sequent cloud expansion. Therefore, Dj^ • D R = and 
the interference pattern is that of Fig. [3] More gener- 
ally, for a given magnetic field gradient one can choose 
the field each cloud is subjected to by adjusting the rel- 
ative distance between the clouds, and control the pair 
wavefunctions (or D vectors). 



V. ANISOTROPIC EXPANSION 

Until now, we considered Fermi condensates trapped 
in a harmonic potential only in the strongly interacting 
BEC limit, where fermions form tightly bound molecules 
and the internal degrees of freedom of the fermion pairs 
do not play an important role. However, if one moves 
away from the BEC limit towards unitarity, the average 
pair size increases with decreasing interaction strength, 
and the internal structure of fermion pairs can dramati- 
cally change the condensate properties when the pair size 
becomes comparable to the inter-molecular spacing. 

In this section, we discuss the expansion of a harmon- 
ically trapped cloud away from the BEC limit (but still 
on the BEC side of the Feshbach resonance). In such 
case, the method used to derive the effective action Eq. 
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([5]) is not directly applicable, since the internal-degree- 
of-freedom wave function is no longer localized within a 
small volume. However, if the trapping potential C/ ex t 
varies slowly in comparison to the coherence length, a 
semiclassical approximation can be applied. Within this 
approximation, one can first derive the effective action 
in a free space (without trapping potential) , and add the 
potential afterwards. Again, we consider here only two 
cases: a non-unitary case where fermions are in a single 
hyperfine state, and unitary cases where fermions are in 
two equally populated hyperfine states. 

In such cases, the quadratic term of the effective action 
takes the form 



? (2) 
3 cff 



JdtJ dRDj^R,*) -£D m ,(R, t), 



(20) 



" E y -ctV,V J /(4m) + 2C/ ext (R) + 
we use the basis of spherical 
and define Dj me by dj(r,r') = 



where C = a me 
id mt d t . Here, 
harmonics Y\, mi 

Em 4 Dj jme (R)r) 0t i(p)Yi jme (p), and assume that the or- 
der parameter is dominated and characterized by the 
spherical harmonics with either Yj., mf= o or Y\ tmi —±\, 
The coefficients a, , and d in Eq. (f2T)]) can be obtained 
by considering the free space problem and transforming 
it into momentum space [12]. For a weak trapping po- 
tential, the coefficient becomes 



E 



X(k) (3Y(k) 



4E 2 (k) 16E(k) 



where E(k) = f 1)k + £ 2 .k, £ Q ,k = k 2 /2m - fj, a , 
X(k) = tanh(/3£i, k /2) + tanh(/?£ 2 , k /2), and Y(k) = 
sech 2 (/3^i jk /2) + sech 2 (/3^ 2 ,k/2). The symmetry function 
cj)(k) is denned by V(k, k') = / dpV(p) exp[i(k-k')-p] = 
Vcj)(k)<j>(k')Yi, me (k)Yj* m , (k'), and the angular average is 



dkkikjYi me (k)Yi (k) 



(22) 



For the case where mg — 0, this angular average is 



anisotropic with Kq 



yy 



I and Kg 31 = 3/5. Here, 



we choose the x direction to be the quantization axis. 
Therefore, the coefficient c 1 q (which is directly related to 
the Ginzburg-Landau coherence length £„■) is diagonal 
and anisotropic, hence acquiring a mass anisotropy Mj = 
2m /cq. This mass anisotropy reflects the higher angular 
momentum (p-wave) nature of the order parameter for 
paired fermions, and it is absent for s-wave Fermi and 
atomic Bose condensates. 

As an example, we consider a p-wave Fermi condensate 
in an axially symmetric trap where uj x — u) v 3> UJ Z , with 
a magnetic field applied along x (chosen as the quan- 
tization axis) to tune through the Feshbach resonance 
(see Fig. 01 . Since the resonances for different mi states 



(a) y 





FIG. 4: p-wave Fermi condensate of p^-symmetry (me = 0) 
(a) in an axially symmetric trap and (b) upon release from 
trap. Notice that the axial symmetry is lost in the xy plane 
due to the anisotropic effective mass (interaction). 



are split as in 40 K [Tj, [T3], it may be possible to ad- 
just the magnetic field such that fermions are paired in 
the mi — (jpx) state only. In this case, the p-wave in- 
teraction leads to the formation of p x symmetry pairs, 
which are more strongly correlated along the x direction 
(6r > iy = & or M x < M y = M z ). Thus, it is eas- 
ier to accelerate the cloud along the direction of lighter 
mass M x such that the cloud expands faster along the 
x-direction than along the y-direction, hence breaking 
the axial symmetry. This anisotropic expansion due to 
p-wave interactions also occurs for a completely isotropic 
trap, and it is very different from the anisotropy inver- 
sion (in the xz and yz planes) found in axially symmetric 
traps for s-wave Fermi condensates [331 ] . The anisotropy 
inversion is related only to the anisotropy of trapping po- 
tential, while the anisotropic expansion discussed here is 
due to anisotropic interactions. 

In Fig. [5j we show the cloud anisotropy ratio tl = 
L x /L y as a function of the effective mass anisotropy ratio 
tm = M y /M x = £ 2 /£y. The anisotropy effect disappears 
in the BEC limit as the effective masses become isotropic, 
but becomes more evident towards unitarity. The values 
of tm change as a function of the scattering volume a p , 
and vary from tm = 1 in the BEC limit (a v — > + ), to 
rjj = 3 in the BCS limit (a p — > 0~). Since our theory is 
valid only on the BEC side where the fermion chemical 
potential p < 0, the maximal theoretical anisotropy is 
reached near p = (which is also close to the unitar- 
ity limit a p — > ±oo [l^]), leading to a 10% anisotropy 
(tm ~ 1.1) for trapped 40 K in the p^-state (mi = 0). 
Investigations on the BCS side (fi > 0) and at unitar- 
ity require the inclusion of Landau damping which leads 
to the decay of Cooper pairs, and are beyond the scope 
of the present theory [341 ]. However, even in this regime 
fermion pairs survive at least initially in time-of-flight 
experiments for s-wave systems [35| . 
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FIG. 5: Cloud anisotropy ratio tl = L x /L y as a function 
of effective mass anisotropy ratio tm = M x /M y at time A 
(Solid lines). Dashed line indicates the saturated behavior at 
A — > oo. 



VI. CONCLUSIONS 

In summary, we considered Fermi condensates with s- 
and p-wave interactions (e.g. 6 Li and 40 K), and derived 
the equation of motion for the p-wave case in a vector 
boson representation on the BEC side of the Feshbach 
resonance. We derived general equations of motion for 
the order parameter in the unitary and non-unitary p- 
wave cases, and showed that the equation of motion can 
be simplified to a similar form as time dependent Gross- 
Pitaevskii (TDGP) equation for vector atomic Bose sys- 
tems, when the superfluid consists of fermions in a single 
hyperfine state, or in two equally populated hyperfine 
states. Within these two special cases, we described the 
time evolution of the vector order parameter approxi- 
mately by scaling the spatial coordinates, and found that 
p-wave Fermi condensates behave very different from s- 
wave Fermi condensates in the following aspects. 

First, the matter- wave interference and the Josephson 
effect of two p-wave Fermi condensates has an angular 
effect due to the vector nature of the order parameter. 
When the dot product of the vector order parameters of 



the left and right condensates reaches its maximum, the 
interference pattern is most visible and Josephson current 
is the largest. However, when the dot product of vector 
order parameters is zero, the interference pattern disap- 
pears and the Josephson current vanishes. This effect is 
absent in the BEC limit of s-wave Fermi superfluids, as 
well as in scalar Bose systems. It was also proposed that 
the relative orientation of the order parameters of two 
p-wave condensates can be controlled by applying a con- 
stant plus a gradient magnetic field which makes clouds 
cross different Feshbach resonances and prepares the vec- 
tor order parameters in orthogonal or non-orthogonal 
configurations. These general considerations were then 
applied to the specific case of 6 Li, where Feshbach reso- 
nances are split into |11), |12) + |21), and |22), depending 
on the different hyperfine states of the colliding atoms. 

Second, we showed that anisotropic p-wave interac- 
tions lead to anisotropic effective masses for a given 
orbital symmetry, as unitary is approached from the 
BEC regime. Furthermore, for cigar-shaped clouds with 
axial symmetry, we found that the cloud expansion 
is anisotropic in the radial plane and expands more 
rapidly along the direction of smaller effective mass due 
to the anisotropy of p-wave interactions. We empha- 
sized that this anisotropic expansion is a result of the 
anisotropic p-wave interactions and occurs even in a com- 
pletely isotropic trap, in sharp contrast with the standard 
anisotropy inversion observed during expansion due to 
anisotropic cloud confinements. In addition, we would 
like to stress that the anisotropic expansion should oc- 
cur not only for p-wave, but also for any higher an- 
gular momenta (d-wave, /-wave, etc..) Feshbach reso- 
nances. Lastly, we pointed out that the orbital symme- 
try of the order parameter for p-wave condensates can 
be directly probed through cloud expansions and that a 
potential candidate for such experiments is 40 K, where 
p-wave Feshbach resonances are split depending on the 
internal angular momentum states (£ = I, me = 0) and 
(£ = l,m e = ±1). 

We would like to thank NSF for support (Grant No. 
DMR-0304380). 
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